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We study the geometry of the event horizon of a spacetime in which a small compact object 
plunges into a large Schwarzschild black hole. We first use the Regge- Wheeler and Zerilli for- 
malisms to calculate the metric perturbations induced by this small compact object, then find the 
new event horizon by propagating null geodesies near the unperturbed horizon. A caustic is shown 
to exist before the merger. Focusing on the geometry near the caustic, we show that it is deter- 
mined predominantly by large-i perturbations, which in turn have simple asymptotic forms near the 
point at which the particle plunges into the horizon. It is therefore possible to obtain an analytic 
characterization of the geometry that is independent of the details of the plunge. We compute the 
invariant length of the caustic. We further show that among the leading-order horizon area increase, 
half arises from generators that enter the horizon through the caustic, and the rest arises from area 
increase near the caustic, induced by the gravitational field of the compact object. 
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I. INTRODUCTION 

The forthcoming prospects for gravitational-wave de- 
tection [T] have motivated a great deal of study into black 
hole mergers, using both perturbative |3] and fully 
numerical methods [4]. The most immediately relevant 
results are the gravitational waveforms emitted by such 
mergers, since accurate theoretical templates for these 
waveforms are crucial for their first detection [5]. As a 
further step, however, studying features of spacetime ge- 
ometry in the strong-field region, as well as their possible 
connection with out-going gravitational waves, will facil- 
itate the use of gravitational waves as a tool for studying 
the non-linear and highly dynamical regime of geometro- 
dynamics [BHSj- 

In this work, we focus on the defining geometrical fea- 
ture of a black hole spacetime, the event horizon — 
a "surface of last return" which separates those points 
which can be connected to future infinity from those 
which cannot jlOl I16j . The horizon is a three-dimensional 
surface in a four-dimensional spacetime, it is a null 
surface, traced out by a two-parameter family of null 
geodesies. At each point V on the event horizon, its tan- 
gent space T-p contains at least one null generator, which 
is tangent to a null geodesic whose entire future devel- 
opment lies on the horizon. There also may exist points 
on the horizon, e.g., Q in the figure, at which two or 
more null vectors lie tangent to the horizon. Although 
future developments of these null generators will stay on 
the horizon, their past developments do not belong to the 
event horizon. Points like Q are called caustic points. 

By choosing a time slicing, one can take a three- 
dimensional cross-section of the spacetime to get the hori- 
zon's geometry at the present time ~ a two-dimensional 
surface in a three-dimensional space. To state an exam- 
ple, the event horizon of a static Schwarzschild black hole 
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FIG. 1: Spacetime diagram of the merger of two black holes. 
The horizontal cross sections give the three-dimensional ge- 
ometry of the event horizons as time progresses. The null 
rays which trace out the horizon are given by the black and 
red directed lines. The black lines originate from the horizon 
at past infinity, while the red lines enter the horizon through 
the caustic. A caustic point Q is shown with its two null 
generators entering the horizon. 



is a cylinder (S^ x R) in spacetime, or viewed in terms 
of its time slicings, it is a spherical surface (S"^) in space. 
The horizon of a black-hole merger spacetime, on the 
other hand, has a "pant-leg" shape, which, in terms of 
its time slicings, looks like two roughly spherical objects. 
These are the horizons of the merging black holes which, 
evolving as a function of time, merge to form a larger 
spherical object. These time slicings of the horizon are 
often referred to as the event horizon itself. A general 
operational way to obtain the event horizon, e.g., for the 
binary black-hole merger spacetime, is to first go to the 
final state of the spacetime, in which a final, nearly qui- 
escent, black hole exists, with an easily identifiable, late 
portion of the event horizon. Null rays on that horizon 
can be propagated backwards in time, and trace out the 
entire event horizon P^15j . 



As it turns out, the existence of a caustic on the 
event horizon is quite general, most notably in black hole 
merger spacetimes. Intuitively, a caustic develops before 
the merger because null generators, propagating back- 
wards in time along the horizon of one of the merging 
black holes, are gravitationally Icnsed by the field of its 
companion, causing some of them to cross and necessar- 
ily leave the horizon, as is shown in Figure [Tl For black 
hole merger spacetimes, numerical simulations confirm 
this intuition [11] [121 US] , and moreover the mathematics 
of general relativity requires that a caustic form under 
generic merger conditions |17) . 

Although perturbation theory has been applied to 
study the deformation of the event horizon due to tidal 
fields dHi 19, 22 24], as well as dissipation caused by this 
deformation [20, i2lj [^5H?7] , these previous works did not 
consider caustics on the horizon caused by plunging ob- 
jects. In this paper, we study the caustic due to the 
plunge of a small point mass (with mass-energy fi) into 
a big non-spinning black hole (with mass M), using per- 
turbation theory. 

After working out the metric perturbations, we 
locate the new event horizon by propagating null 
geodesies around the future horizon of the unperturbed 
Schwarzschild spacetime. Since the perturbation we ap- 
ply is only valid to within a distance S> /i from the 
small black hole (e.g., as measured in its local asymp- 
totic rest frame), as we shall see later in the paper, our 
study will have to exclude a region in the event hori- 
zon that has area of approximately 0(/i^) towards the 
final future horizon. Fortunately, most of the change 
in geometry is caused by rays that travel at a distance 
^ y/fiM ^ ^ from the small hole, where gravity is still 
weak. More specifically, since the final hole will have a 
mass of M[l-|-/i/M -|-0(/i^/M^)] and angular momentum 
of 0(^M), its area wiU be 16n^AP[l+2iJ,/M+0{^i^/AP)] 
— linear black-hole perturbation will account for the 
leading order area increase due to the plunge of the black 
hole, in particular the contribution due to those rays that 
enter through the caustic. Moreover, as it turns out, rays 
significantly influenced by the small black hole all tend 
to go close enough to the small black hole, such that 
the small black hole's influence can be approximated as 
"instantaneous". This allows us to develop an impulse 
approximation (also called the Born Approximation in 
scattering theory 155") that leads to an analytic descrip- 
tion of the caustic and geometry around it. 

This paper will be organized as follows: In Sec. [ll| 
we first briefly review the Regge- Wheeler and Zerilli for- 
malisms for black-hole metric perturbations, and then 
apply these results to the propagation of null geodesies 
near the future horizon, writing down their evolution 
equations, which are driven by "forcing terms", which 
can in turn be written in terms of Zerilli and Regge- 



Wheeler functions. In Sec. Ill we develop an impulse 
approximation for null-ray propagation and show that 
the deflection of null rays depends only on the time inte- 
gral of the forcing terms, instead of their detailed de- 



pendence on time. In Sec. IV we apply the impulse 
approximation to the propagation of null rays close to 
the unperturbed horizon, obtaining the caustic structure 
of the new horizon. In Sec. [V| we calculate the event- 
horizon area change due to new rays that enter the hori- 
zon through the caustic and due to increase in area in- 
duced by the small object. In Sec. |VI[ we relax the point- 
particle assumption and treat the case of one-dimensional 
"strings" falling into the black hole. This results in sig- 
nificant changes to the caustic geometry and the distri- 
bution for area increase. In Sec. |VII[ we summarize our 
main conclusions. 



II. METRIC PERTURBATION AND 
DEFORMATION OF THE EVENT HORIZON 

A. Metric Perturbations: the Regge- Wheeler 
Gauge 

Suppose /i <^ M, and using the geometrical units of 
G — c = M — 1, we consider a background Schwarzschild 
spacetime with the black hole mass set to unity: 



ds^ = -(1 - 2/r)dt^ + (1 - 2/7 



'dr' 



^dn^ 



(1) 



The small black hole's world line in this background 
spacetime, up to leading order in /i, is a timelike geodesic. 
We consider a first-order perturbation (at 0{p) order) in- 
duced on this background spacetime by the small black 
hole. Dependence of the metric perturbations on angu- 
lar coordinates and (j) can be decomposed into scalar, 
vector and tensor harmonics, and classified according to 
parity. As shown by Regge and Wheeler [2^ , a choice of 
gauge (the Regge- Wheeler gauge) allows us to eliminate 
all but 6 fields, and write 



dsl 



[(1 - 2/r)i/rdi 



2m 



dtdr+ 
iir(l - 2/r)-'dr^ + r^K^'^dn^] F'" 
-K2/i[,™ [smeY^g"'dtd<j>~ CSC OY^^dtde] 
+2h["^ [sin 6Y^^drd(j) ~ esc dY^^drdO] (2) 

Here {Hq,Hi,H2,K) are "even-parity" perturbations 
with a parity of (—1)', and {ho, hi) are "odd-parity" 
perturbations, with parity of (—1)'+^. Henceforth in 
the paper, we drop the (Im) dependence of all metric- 
perturbation fields. 

Regge and Wheeler |29] deduced 10 linearized Einstein 
Equations for these 6 fields. Among these, 7 are even- 
parity and 3 are odd-parity. Zerilli [31^ showed that the 
monopole and dipole perturbations can be found exactly 
by integrating the equations of motion. For multipoles 
with I > 2, however, it is not clear from the Einstein 
equations whether one can solve for the metric pertur- 
bations systematically, e.g., as an initial-boundary- value 
problem. However, for odd and even-parity perturba- 
tions respectively, Regge-Wheeler [21] and Zerilli [30] 
were able to construct functions Q (often referred to as 



the Regge- Wheeler function) and Z (often referred to as 
the ZeriUi function) that satisfy wave equations in vac- 
uum, which can be solved using standard approaches — 
and all metric perturbation fields can then be expressed 
in terms of Z and Q. Zcrilli 31 further worked out 
the source terms that appear on the right-hand side of 
the wave equations, when a point particle falls along a 
geodesic — as well as modifications that must be made 
in the relations between Z and Q and metric perturba- 
tion fields. By these procedures, we can find the metric 
perturbations for all I values. 



B. Perturbations with I > 2 

We briefly review the procedures used to obtain these 
perturbations in Appendices |A] (for even parity) and |B] 
(for odd parity) . Here we simply list the conversion equa- 
tions: 



We have also defined 



n =r + 21n(r/2-l). 



(14) 



which is often referred to as the tortoise coordinate. 

Since we are interested in the deformation of the fu- 
ture event horizon, we specialize the relations (3]l8) to 
the event horizon, where (r — 2) — > and Z and Q are 
functions oi v = t + r.. alone: 
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where 
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and the wave equations, 

d^z d^z 

dri dt"^ 
(PQ__(PQ_ 
dri df^ 



oZ 



-Vi^{r)Z 
-V^{r)Q = S\ 



(9) 

(10) 
(11) 



The quantities Shq through SV, are placeholders for the 
source terms in Equations ( |A11[ ), ( |A21|A23[ ), and ( [BS" 



B9|). The source terms S^^ and S^^^ take rather cumber- 



some forms and have been consigned to Equations ( A20 ) 



and (B7| of the Appendix. The potential terms in the 



ZeriUi and Regge- Wheeler wave equations are 

2\ A2(A + l)r3 + 3A2r2 + 9Ar + 9 



Vf{r)^2[l- 



r3(Ar + 3)2 



V^{r)={l- 



2(A-|-1) 6 



(12) 
(13) 



where Q{v) is the Heaviside step function, ^ = rrioE is 
the mass-energy of the small hole in the Schwarzschild 
background (toq is rest mass and E the specific energy 
associated with its geodesic world line), and YJ^ = 
yj^(0Oj 4'o) is the complex conjugate of the spherical har- 
monic taken at angular coordinates (^Oi^o) where the 
particle meets the horizon. 

As we can see from Eqs. ( 15p8| , the metric pertur- 
bations are given by ingoing waves (i.e. functions of v) 
multiplied by various powers of (r — 2). We may write 
this dependence more explicitly as: 



H^ = 



K = IC{v) Ho = 
Hiiv 



H.^ 



r-2 

ho = f]o(«) hi = 



Hojv) 
r-2 

H2[v) 
r-2 

M^ 

r-2 



(19) 



where K, Ho, Hi, H2, i}o, and f)i are defined to make 
( [191 ) consistent with pSflSl ). 



C. Low Multipole (I < 2) Perturbations 

For low values of I, not all Regge- Wheeler fields are 
involved in parametrizing the full metric perturbation. 
The linearized Einstein Equation, consequently, will be 
dramatically simplified. For example, when / = 0, vector 
and tensor harmonics all vanish, while for I — 1, only the 
tensor harmonics vanish. A direct mathematical conse- 
quence is that a wave equation cannot be constructed for 
perturbations with I < 2 — yet equations are simpler 
so that they can be solved directly. In Appendix G of 
Ref. [31], Zerilli provides solutions to all these cases for 
a point particle perturbing a Schwarzschild space-time. 
We do not repeat his derivation, but merely state the 
results. 



1. Monopole {I = 0) Term 

The monopole perturbation [I = 0) is associated with 
the mass of the black hole. At a distance r, it is related 
to the amount of mass enclosed within the coordinate 
sphere with radius r. Since vector and tensor harmonics 
do not exist for I = 0, the only perturbation fields are 
Hq, Hi and H2- As Zerilli has shown in Appendix G of 
Ref. [3T] , after a gauge transformation, the only surviving 
even-parity terms are Hq and H2 ■ For a plunging particle 
of mass-energy /i, we have: 



Hn — Ho — 



2V47r/i 



e{v) 



or equivalently, 



(20) 



(21) 



•Ho = -Ha = 2V4^^e(t;) 
where Q{v) is the Heaviside step function. 

2. Odd-Parity Dipole {I = 1) Term 



The odd-parity dipole perturbation represents the spin 
imparted by the small black hole's orbital angular mo- 
mentum. This is slightly less trivial than the / = case, 
but can be simplified by a gauge choice which makes hi 
vanish and /iq approach a constant in time both before 
and after the plunge. A gauge can be chosen in such a 
way that metric perturbation before the plunge vanishes, 
and after the plunge acquires a value that depends on the 
orbital angular momentum L of the plunging particle: 



5gtr -- 
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(22) 
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2 r . . , 


- Ly cos 0] 
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[—Lz sin^ 6 + {Lx cos 6 
'Ly sin 9) sin 6 cos 6] 



(24) 



This corresponds to a metric perturbation of: 
f)o'" 



. Im 
"0 



y^^?!^!2^e(i-T(r)) 



v/47r/3moL"e(i-T(r)) 



(25) 



where L™ is the spherical-harmonic representation of the 
plunging orbital angular momentum L - i.e. L'^ = U' and 

i±i=2-i/2(L,±zL,) 



3. Even-Parity Dipole Term 

In the absence of a source, the even-parity dipole term 
is a gauge that can be eliminated by transforming to the 
center-of-mass frame, in which the center of mass lies 



at rest at the origin of the coordinate system. When 
a source is present, however, the term cannot be elim- 
inated, since it is not possible to gauge away a source 
term (although a gauge tranformation can concentrate 
the even-parity term along the path of the particle). 
Choosing to work in a gauge where i^iT = 0, the Einstein 
Equations can be integrated exactly to give: 

/i™(0+'^'/(™Wq(^_^(^)) (26) 



Tjlrn 



Hi"" = 
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(r-2)2 



(r-2)2 

flmv) 



Q{r-R{t)) 
Q{r-R{t)) 



where 



h„,{t) = ^TTli{R{t)-2)Y, 



2) 



Im 



(27) 
(28) 

(29) 



-t/2 



SO /{^ 



Near the horizon, we have {R{t) 

-\fim and f'{^ = \fini- In Section (|IID|), where we 



trace out the structure of the perturbed horizon, we will 
make use of the values of the metric coefficients, as func- 
tions of V, along lines of constant u = t — r^ (where 
(r - 2) - e"/4 and {R{v) - 2) ~ e""/"*). Doing so here, 
and setting i; = to correspond to the point at which 
R(v) crosses the line of constant u (so r — 2 — R{v) — 2 
at w = 0), we find: 

njlm njlm njlm 

rtg — til — '^2 



= %'K^ie-''/'^Q{v)Y*, 



(30) 



Of importance, we note that all metric coefficients vanish 
at past and future infinity, indicating that, in these limits, 
the coordinate frame is centered around the large black 
hole. Since it is the large black hole that we are interested 
in, this is the proper coordinate frame to use. 



D. The Deformed Event Horizon 

In order to analyze the horizon deformation caused by 
metric perturbation fields ( 15p8 l, we need to study the 
propagation of light rays near the horizon. We will do so 
in the light-cone Kruskal-Szekres coordinates, which offer 
the distinct advantage of having non-singular light cones 
around the horizon. The Kruskal coordinates {V,U,d,4>) 
are related to their Schwarzschild counterparts by: 



where 



X" = V = e-/* 


(31) 


X^ = U = -e-"/4 


(32) 


x^ = e 


(33) 


X^ = cf) 


(34) 


V = t + r.j. 


(35) 


u — t ~ r^, 


(36) 



We model the unperturbed horizon as a set of null gener- 
ators, parameterized by V. To distinguish these horizon 
generators from the rest of the null rays in the system, 
we impose the final condition [/ — >■ as ^ — >■ oo. The 
generators of the unperturbed horizon are then given by: 



V{V) 
U{V) 
9{V) 
0(T/) 



X%V) = V 
X\V) = 
X^{V) = const 
X^{V) = const 



(37) 
(38) 
(39) 
(40) 



These generators satisfy the geodesic equation, modified 
to account for the non-affine parametrization: 



d^X'' _ -^ dX" dXP 

dV^ ^^ ''p~dy~dv 



9- 



dXt" 
W 



(41) 



Here, V^ refers to the Kruskal Christoffel symbol. The 
non-perturbed horizon at [/ = indicates 



.9 = rSJo = 



rf,o = o. 



(42) 



Note that the unperturbed event-horizon is affine 
parametrized. 

On a perturbed metric, the rays themselves will be 
perturbed; supposing we still parametrize the horizon by 
y, then we need to modify X^ -^ X' + SX^ (j = 1, 2, 3) 
and g -^ g + Sg. To first order in the metric perturbation, 
we obtain the following equations of motion: 



Sg - ~2f 



d{SX^) 
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dV^ 



-2T 
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(43) 

d{5X') 

dV 

(44) 



are null rays, and always will be no matter what space- 
time they propagate through. Given 9{V) and (t>{V) from 
(45), we can derive an equation for 5U{V) by setting 



5^1' "f 



= 0: 



1 d{5U) 
e dV 



1 
16' 



goo 



m 



sm 



fdm\ 



(46) 

The equations in their present form are solvable but cum- 
bersome. One can rewrite them in a more intuitive form 
by transforming into the ingoing Eddington-Finkelstein 
form, substituting U for r and writing V in terms of v. 
Upon simplification, equations (45M6I transform into: 



^4)---- 


6^ 


+ 4>' 


. 9 

sm 
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dv AJ dv 








(48) 


dv AJ dv 










(49) 



where A means dA/dv for any A, and F^ , F^ , and F'^ 
are "forcing terms" that arise from the small black hole's 
perturbing field. 



E. The Forcing Terms 

The "forcing terms" F^, F^, and F'^ tell the rays on the 
horizon how far to bend in the small hole's gravitational 
field. In terms of the Kruskal-frame metric perturbations 
and Christoffel symbols, they are given by: 



Note in this case that since Sg ^ 0, the perturbed event 
horizon is no longer affine-parameterized by V^. A careful 
derivation shows that all of the Christoffel symbols in ( 42 ■ 
44 1 are finite in the limit r — )■ 2. This is necessary in order 
for the perturbation theory to be well-posed. Moreover, 
most of these Christoffel symbols (but not their pertur- 
bations) vanish outright, yielding the following equations 
of motion for the angular coordinates 0{V) and 0(T^): 



d^Sd) 
dV^ 



°^ 00' 



d^iS(t>) 
dV^ 



"^ 00 



(45) 



In principle, we could have also derived an equation 
for SU{V). However, as we will see in Section |IV[ in 
the region of interest, i.e. near the caustic, the angular 
perturbations will scale as 0(/^^/^), while the radial per- 
turbations scale with the higher power 0{fi). It will not 
suffice to compute the geodesic equation to first order in 
the perturbation; in a rigorous treatment, higher-order 
perturbation terms would be needed here. 

Fortunately, there is another way to derive the radial 
equation that does not involve a cumbersome higher- 
order perturbation expansion. The horizon generators 
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Sg 



-^ogoo 

~Sf^ 
16^ °° 

1/2 



16 



-ST 



00 



(50) 
(51) 
(52) 



While the Christoffel symbol is not a tensor, its pertur- 
bation (5r[J is. The Kruskal-coordinate components of 
this tensor are related to its Schwarzschild components 
by the coordinate transformation: 



8X>^ (9t^ Br'^ 

srt" = sr" 
''p ax" '^'-'dx'^dxp 



The Schwarzschild Christoffel symbols 
1 



(53) 



STL 



9'^"{S9<yv-p + Sgap-v - Sg^p-a) (54) 



are calculated using the metric i n (pl a nd substituting 
the horizon metric perturbations ( 15||18). This leads to 



the following forcing terms: 
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(55) 
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(59) 



In terms of the Zerilli and Regge- Wheeler functions, we 
may write this as: 
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III. IMPULSE APPROXIMATION 

Up to this stage, a straightforward calculation using 
the techniques developed above can already compute the 
metric perturbations, the shape of the horizon, and hence 
the structure of the caustic. This involves numerically 
solving for the Zerilli and Regge- Wheeler functions for 
particular geodesies, and inserting them into Eqs. (47- 
49) and ( 55][61 ). In particular, the wave equations can 
be solved either in the time domain, or in the frequency 
domain. For example, for Z, we have 






-V,^ 



Z{uj,r) ^Sf^, 



with out-going boundary condition 

Z(a;,r,)~e±''^''*, r, -> ±00 . 



(62) 



(63) 



However, when we look on the short lengthscales that 
characterize the caustic, or conversely if we look on the 
long timescales that characterize the initial and final 
states of the big black hole, the results can be greatly 
simplified by approximating each forcing term as an in- 
stantaneous impulse. 



A. The Approximation 

Consider a generator on the large black hole's hori- 
zon which is deflected by the small hole's gravitational 



field. Let 9 denote the impact parameter of this scatter- 
ing problem. Most of the deflection will occur within a 
time interval 

T ^9 (64) 

Define a "boundary" 9), such that ,JJi ^ 0f, <C 1. All 
rays with 9 < 9b are considered to be in the vicinity of 
the caustic, while rays with 9 > 9^ are considered to be 
in the hulk of the horizon, 9b functioning as the bound- 
ary between these two regions. We consider the impulse 
approximation for both cases separately: 



1 . Vicinity of ttie Caustic 

As we will show in the next section, the geometry of the 
caustic depends only on those rays with 9 < yjfi <^ 9b - 
that is, rays in the vicinity of the caustic. The deflection 
timescale for these rays is thus ^ 1, and the deflection 
can be approximated as an instantaneous impulse. This 
is the intuitive essence of the impulse approximation; we 
discuss it in more rigorous detail below. 

we show that Z and Q 



In Appendices |C 1| and C 2 



satisfy scaling relations for large I, namely: 

The forcing terms (601161) depend on the Zerilli and 
Regge- Wheeler functions, and therefore satisfy similar 
scaling relations. As an illustration, the scaling rela- 
tion for the even term f^^ is plotted in Figure bl As 
I increases, the forcing term increasingly resembles an 
impulse. The impulse approximation thus replaces /j 
and /, ° with delta functions in time: 
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f1:!S{v) (67) 



2. The Bulk 



We will prove in Section IV that 



^ ^/9: let us for 

The deflection of 

69 <^ 9. Taking 



now assume a priori that this holds, 
rays in the bulk is then infinitesimal 
this to be the case, the forcing terms in ( 47p9 ) lose their 
dependence on the deflection and become functions of 
time alone. We can then integrate ( 47p9 ) to give: 
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values 2 < Z < 24 

denoted by darker lines. The red dashed line is the empirical 

limit curve Svre"^'^; the significance of this curve is touched 

on in Section Till CI 



2. Odd Parity, I > 2 



f(o) 



Here, we integrate /;^ over time with the help of (|61|), 
to obtain: 



j(o) 
•f Im 



Q\ 



+CX3 



Q{v')dv' 



(74) 



Unlike the Zerilli function, the Regge- Wheeler function 
approaches zero as w — >■ ±cx), so the Q term in ( 74 1 drops 



out. The time integral of Q, however, does not vanish, 
and as we show in Appendix |C2[ this term scales as 



1-^y: 



m 



Thus, for large I, 



j(o) 
•' Irii 



r^Y, 



.(0) 



im,e 



r^Y, 



'(0) 



Im 



(75) 



We need not compute the precise form. What matters is 
that the odd term (75) scales as l~^ , while the even term 
scales as l~^ . This means that for large values of Z, the 
odd terms may be neglected and only the even impulse 
need be considered. 



On long timescales \v\ ^ 1, the angular perturbations 
depend only on the time integral of the forcing term; be- 
cause of this, we can replace the forcing term with an 
instantaneous impulse. The radial term, conversely, dep- 
nds on the value of F^ &i v = ±oo; as we will show, this 
gives rise to an increase in the black hole's radius by an 
amount 5r = 2/i. 



B. The Impulse 



Even Parity, I > 2 



Integrating /^^ over time with the help of (|60|), we 



find: 



-(e) ^ 32^/zy;^(°) 

■''"* 2A + 3 



dv 



+00 



(71) 



The Zerilli function approaches a constant as u — >■ ±oo, so 
the dZ/dv term drops out, but the jump in Z contributes 
to the final result. In Appendix |C ![ we show that the 
jump is given by: 



AZ 



l+OO 



SiTfirY* 



(A + l)(Ar + 3) 



(72) 



Substituting this into (71 1, we find an impulse of: 



-(e) ^ WnnYi 

J Im 



.(0) 



Irn 



A + l 



Note that this scales as l^^Y^^^^ for large /. 



(73) 



3. Monopole and Dipole Terms 

The monopole term accounts for the mass increase of 
the large black hole, which, to first order, is /i. By the No- 
Hair Theorem, we infer that on long timescales \v\ ^ 1, 
it gives rise to the following radial perturbation: 



Sr 



2/1 (u — )■ oo) 
(w -> 0) 



(76) 



Since the monopole term is isotropic, it does not give rise 
to angular perturbations. 

The odd-parity dipole term likewise cannot be viewed 
as an instantaneous impulse, but instead should be 
treated as a constant forcing term on long timescales: 



/i 



tl = ~lV^7^moL"^e{v) 



(77) 



Using Eqs. ( 49 ) and (571, we see that this results in a slow 



rotation at a rate (j){v) — gL (if L is chosen to point along 
the z-axis). Recalling that the horizon perturbation is 
done along lines of constant u (for which dt = 2dv), we 
see that, in the Schwarzschild frame, the rotation rate is 
nothing more than (Ikcii = L/A for slowly rotating Kerr 
black holes. 

Unlike its odd counterpart, the even-parity dipole term 
may be treated as an impulse. Substituting (30) into 
(58), it takes the following form: 



ftl = 8nt,e-^/'Q{v)Y,: 



(0) 



(78) 



On long time scales, this resembles an impulse of the 
form: 



m - l6n^,YZ" = 



(0) _ IGTTfiY^. 



*(o) 



A + l 



(79) 



This agrees with Equation ( 73 1 , which we calculated only 



for I > 2. Thus, (73) is valid for all / > 1. 



C. Shape of the Impulse 

In this section, we have shown that the forcing terms 
can be approximated by an instantaneous impulse. Near 
the caustic, this approximation is certainly valid, as we 
show in Appendix IdI so the shape of the impulse has no 
first-order effect on the caustics; however, it will prove 
enlightening to consider its shape nonetheless. 

As we showed in Figure |2] as we increase the value of 
I, the forcing term approaches the following scaled limit 
curve: 



/(:U8.Mr^K:re-'i^i/^ 



(80) 



In previous sections in this paper, we approximated this 
as a delta-function. Here, we choose to retain the time 
dependence and proceed directly to calculate the effect 
on the horizon generators. Without loss of generality, as- 
sume that the small black hole falls into the horizon at 
9 = 0. Then only terms with m = need be considered. 
Restricting ourselves to the region <C 1, we need con- 
sider only those terms with large I values. The forcing 
terms at {9,(f)) will be given by: 






16 961 ^^-^ 



lo\v)Yio 



1 1 



16 sin^ ' 



I 



f^o\^)Yio = 



(81) 
(82) 



Proceeding, we simplify F^ by noting that Y/o rnay be 
written as a Bessel function for <C 1: 



" 2 d9^ '0 



Y, 



io\ 



^l\v\/2 



Y,i4iie)e- 



l\v\/2 



(83) 



We recall that the forcing term is very nearly an impulse, 
and therefore F^ will only be significantly nonzero when 
\v\ ^ 1. Taking \v\ <^ 1, we see that the term inside the 
sum behaves fairly smoothly as a function of Z - that is, it 
does not change much if we increase I to l + l, or decrease 
it to Z — 1. Thus, we can replace the discrete sum with 
an integral over I, and evaluate the integral analytically: 

F' ^ f^ r I4{ie)e-'\''\'^dl 



46*2 



e-(kl/2e)?^j.(^)^^ 



(84) 



4 (612 + t,2/4)3/2 

A null generator starting at (6*, (f) can be thought of as 
a light ray scattering off of the small black hole with an 
impact parameter b = 20. The transverse acceleration of 
the ray is given by: 

a^ = 2F'^ 

_ b 

-2a* (^,2 + ^2)3/2 (85) 



Recall that this is the transverse acceleration of null rays 
on the horizon induced by the gravitational field of the 
small mass jjl. If we transform into what is analogous 
to a Fermi normal coordinate system centered around 
the geodesic of the infalling mass, the metric near the 
small black hole will be locally Schwarzschild. Near the 
small black hole, the horizon generators become light rays 
propagating through the (locally) Schwarzschild space- 
time of the normal coordinate frame. For weakly scat- 
tered light rays in a Schwarzschild metric, it has long 
been known that the transverse acceleration is given by 
(85), which is twice that predicted by Newton's theory 



of gravitation. [33] 

As our exercise shows, the perturbation theory gives 
a result ( 85 ) which agrees with the intuitive result we 



would expect if we considered only the local behavior of 
the horizon generators in the vicinity of the small black 
hole, on spatial and time scales much smaller than the 
large black hole's radius of curvature, so the large black 
hole's field does not affect the result. Near the small 
black hole, the impulse is the same as it would be if the 
large black hole had been absent. 



D. Results and Accuracy 



Substituting ( |73[ ) into ( 55][57 ), the impulse approxima- 
tion gives the following forcing terms (valid for I > 2 and 
even-parity Z — 1) 






(86) 



y*(0)y; 



F<^ = T:tiS{v)J2 'T f + (odd terms) (87) 



F* = TTfi5{v)Y ^:^ ^ + (odd terms) 



We have not calculated the odd-parity terms explicitly, 
since as we showed in Section III B 2 in the near-caustic 
impulse approximation {9 ^ 1), their effect is negligible 
compared to the even terms. In the far-from-caustic im- 
pulse approximation [9 ^ y/m, \v\ > 1), the odd terms 
will be comparable to the even terms. 

The I = and odd-parity I = 1 terms are an exception 
to our delta-fucntion rule, since they result in a perma- 
nent change in the black hole's mass and angular momen- 
tum, respectively, and cannot be treated as impulses. 

In AppendixlDJ we show that for rays in the vicinity of 
the caustic - that is, for 9 ~ ^/Jl, the relative error due to 
imposing the impulse approximation is at most of o rder 
Of/i^/'*). This is done by splitting the expressions (86 
88 1 into terms with I < ^~^/^ and terms with I > /Li~^, 
showing that the former can be neglected up to a relative 
error 0{ii^^^), and that the latter can be approximated as 
a delta function up to a relative error of at most 0{fi^^^). 



This is only an upper bound, however, and empirically, 
the relative error appears to be of order 0(/^^/^). 

The radial forcing term, by contrast, is not dominated 
by the large-Z terms. As we show in Section |IV[ F^ is 
of order m log m, while the low-/ perturbations give a 
contribution of 0{m). Therefore, the low-/ terms cannot 
be ignored in the radial case as they could be in the 
angular case. However, from ( 55][57 ), we see that the 
derivatives of F^ can be related to the angular terms 
by V°F'' = — 4F(°^j,j >, which is well approximated as an 
impulse. This identity is exact and holds for all /. If 
we then restrict ourselves to a small patch of the horizon 
(say, the neighborhood of the caustic) , F^ will equal its 
large-/ contribution up to a time-varying but constant 
function C^'^' . The impulse approximation for F^ thus 
holds up to a constant term: 



v = -l. 



-0.75 



F'- = -4T:fj,6{v) J2 



Im 



^ Im -^ 

A + 1 



-C'^''\v) (89) 



The constant C^'^\v) depends on the time coordinate, 
but is independent of the and (f) coordinates, falls off to 
zero as \v\ 3> 1, and as we will show, does not affect the 
properties of the caustic. 

Far from the caustic, the deflection of any ray is very 
small - SO <^ - and therefore the impulse approxima- 
tion will be valid at large times up to a relative error 
0{59/6) ^ 1. However, it is important to recall that far 
from the caustic, the impulse approximation only holds 
on timescales \v\ ^ 1, as described in Section HI A 2 
The actual "impulse" is spread out over a timescale 0{1) 
and therefore does not look like an impulse for timescales 
V ^ 0(1). Also, unlike the caustic-vicinity case, in the 
bulk, odd-parity contributions to the impulse are not neg- 
ligible. 



IV. EVENT HORIZON CAUSTICS 

A remarkable feature or black hole merger spacetimes 
is the presence of caustics - regions where the null gen- 
erators on the event horizon cross each other, giving rise 
to discontinuous kinks in the event horizon's geometry. 
While the event horizon may be sharp and pointed in 
the vicinity of a caustic, the metric remains smooth; 
therefore perturbation theory can be applied to study the 
structure of these caustics. The impulse approximation 
in the above section simplifies matters greatly, allowing 
the structure of the caustic to be determined analytically 
as a function of the mass ratio. 



A. Dynamics of Generators that form Caustics 

Let us define coordinates x° = (0, (j)) on the unit 
sphere. Using this notation, we may recast Equations 
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FIG. 3: Event horizon and caustic of a black hole merger with 
m = 0.15, plotted in an ingoing coordinate system at times 
V = —1, —0.75, —0.5, and —0.25 respectively. Black dots 
indicate null generators on the horizon; red dots indicate rays 
which have yet to enter the horizon via the caustic. The large 
black dot is the infalling black hole. 



( |47p9l ) 



as: 



d_ 
dv 



l\ dx" 
4/ lb 



F»(a;^w) 



(90) 



Noting that in the orthonormal basis of spherical har- 
monics 2(A + 1) = — V^ (the Laplacian being taken on 
the unit sphere) and Y*Jf'' — j dUYi*^fi'''^\U), we may 



simplify Equations ( 87p8 ) and (89) to give: 

F'- ^ -4^/x(-2(V2)-i)5(2)(J7),5(w)-|-C(°)(w) (91) 
F° = ^tiV'{-2{V^)-^)5'^'^\n)5{v) (92) 

from which it is clear that 



pa 



4 



(93) 



Therefore, finding the forcing terms boils down to solv- 
ing Poisson's Equation for a point source fj.S^^'^l^). Not- 
ing that the "potential" 



* 



At: 



In [1 - cos 0] 



(94) 



satisfies V^vji = (5(2) (fi) (at least for / > 2, where the 
Zerilli formalism if valid) , the forcing terms take the fol- 
lowing form: 

87r/i*(5(u) + C(°Hu) 

2^1n[l-cos6']5(i;)-|-C(°'(«) (95) 



F"- 



pa 



-27r/xV''*(5(w) 



1 1 + cos 6* 
2^ sin 6* 



S{vK (96) 
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9(v) Caustic Horizon 
^ \ Event Horizon 




FIG. 4: Null event horizon generators, traced back in time. 
The rays inside the caustic horizon cross each other and exit 
the event horizon, forming caustics. The rays outside the 
caustic horizon always remain on the event horizon. 



r{v) = 
0{v) = 



-(8/iV^+)(e"^'-e"/^)+c(«) 



(103) 
(104) 
(105) 



where 



r_ 


= r+ = 2 


(106) 


K 


= -^^,[\og{e+) + c] 


(107) 


e_ 




(108) 



and civ) is defined so that 



(cg is the 6 unit vector on the unit sphere). 

In this equation, we have neglected the I = and odd- 
parity I = 1 terms, which do not function as impulses. 
That aside, ( 95][96 ) give a general expression for the forc- 
ing terms and may be used both in the caustic vicinity, 
and in the bulk, of the event horizon, under the condi- 
tions spelled out in Section |III| In this section, we are 
interested in the caustic structure and accordingly treat 
the case <C 1, where the I = and I — 1 terms play an 
insignificant role. In this regime, ( 47p9 ) take the follow- 
ing form: 



d_ 
dv 



6r = 4^(log(6i) + C)(5(t)) 



d 
dv 



d 
dv 



+Cl»'(i.) 

->> 

const 



(97) 

(98) 
(99) 



where C = —h log(2) is a constant. 



To illustrate the effect of ( 97|[99 l on the null generators, 
let us trace out the path of a given generator. Since the 
generator is a part of the future event horizon, we "start" 
it at future null infinity (w = cx)) and propagate it "back- 
ward" in time. Let r-|_, 0_|_, and (p^ denote the values of 
r, 6*, and the ray "starts" with at u — > -l-oo. Likewise 
denote r_ , f?_ , and (/)_ as the values the ray "ends" with 
at w — > — oo. Equations ( 97]|99 ) are straightforward to 
solve; for w > (after the plunge), we have: 



r{v) 
B{v) 
</>(«) 



0+ 



-civ) 



(100) 
(101) 
(102) 



c'{v)~-c{v)^C^°\v) 



(109) 



These solutions are visualized in Figure [3) Note that, 
as expected, rays enter the event horizon through the 
caustic at 6* = 0. The 2-dimensional time slices of the 
horizon develop kinks at the caustic, in agreement with 
the results of previous literature. [12] Note further that 
the two constants in this result, C and c{v), do not affect 
the internal properties of the caustic, but merely shift it 
in the r-direction. 

For an infalling point mass, equation ([98]) tells us that 
caustics will always form for rays of small enough 9^. 
Naturally, we are inclined to ask: When traced back to 
v — — oo, what does the set of rays that form caustics 



look like? It is fairly easy to see from ( 108 1 that a ray 
will form a caustic if: 



9+ < 2^ 



1/2 



(110) 



Define the caustic horizon as the set of rays with 9^ = 
6c = 2y}/'^ (see Fig. W|. When traced back in time, the 
rays inside the caustic horizon form caustics and leave 
the event horizon, while rays outside the caustic hori- 
zon remain on the event horizon forever. Equivalently, 
rays inside the caustic horizon originate outside the event 
horizon and enter it through the caustic, while rays out- 
side the caustic horizon originate on the event horizon. 



B. Properties of the Caustic 

Globally, the caustic is a spacelike line which lies on 
the future horizon of the black hole. This line can be 
associated with a function 5r{v), defined by solving Eqs. 
( |103|104[ ) subject to the constraint 6* = 0. We find: 



5rc 



-4/i 



log y^4^^(l^^^"^W4^ + ic 



Prior to the merger, the solution takes the form: 



6"/^ + C{v) 
(111) 

The caustic does not lie along a null generator; it is space- 
like and its length is a well-defined and invariant quantity. 
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FIG. 5; The event horizon of the large black hole as a surface 
in spacetime. The small black hole has a mass m = 0.15. Null 
generators are shown as black lines, and the region spanned 
by the generators entering through the caustic is shaded red. 



Taking the line element along the path of the caustic, we 
find: 




FIG. 6: Viewed as a 2-surface in 3-space, the horizon near 
the caustic looks like a cone of angle n — a. This is shown 
for a small cone angle (top) and a larger cone angle (bottom) , 
both cones propagating to the right as time progresses (three 
time slices of the horizon are shown in the figures). Like a 
Cherenkov cone, the speed Vc > c increases as the cone angle 
decreases. 



ds'^ = [-(1 - 2 /re) + Sgyy] dv^ + 2dv dr^ 

pV/2 



f^ 



1 



pV4 



+ 2C^°\v)-2F' 



dv^ 



aV/2 



M 



1 



=11/4 



dv' 



(112) 



Once again, we find that the measurable quantities of the 
caustic do not depend on the "constants" C or C''°-'(w). 



The result (112 1 integrates to a total length 



pi'/4 



%/l - e"/4 



dv — 8y/Jl 



(113) 



Even though the caustic stretches back to v = — oo, it 
has a finite invariant length. 

Locally - i.e. in a local Lorentz frame which includes 
the caustic - the horizon around the caustic resembles 
a cone of angle tt — a propagating outward at a super- 
luminal velocity Vc, as illustrated in Figure [6J At least 
two null rays (fci and ^2, whose spatial components are 
shown in the figure; see also Figure [5]) tangent to the 
horizon meet at the caustic. As we showed previously, 
our caustic has a conical shape, so an entire cone of null 
rays meet at the caustic. To keep things simple, though, 
we only consider two such rays, chosen to be tangent to 
opposite sides of the cone. 

From geometric considerations, we can relate the cone 
angle to the speed of the caustic. We find: 



14 



dx 
~di 



sec{a/2) 



(114) 



where x and t are spatial and time coordinates in the 
local Lorentz frame. Thus, a very sharp cone propagates 
very quickly, while a blunter cone will propagate slower, 
but still superluminally. The speed Vc is in turn related 
to the invariant distance traced out per unit time: 



ds 
di 



^V^ - 1 = tan(a/2) 



(115) 



A third local property of the caustic is the deviation be- 
tween generators at the caustic. If, in a given time slicing, 
we normalize the horizon generators so that fc* — 1, we 
can define a quantity ipc as follows: 



tpc = ^max[|fci -fc2| 

Z ki,k2 



(116) 



The quantity is maximized over all generators (A:i,fc2) 
tangent at the caustic, the maximum being obtained 
when the vectors point along opposite sides of the cone. 
Again, elementary geometric considerations relate ipc to 
the cone angle as follows: 



Ipc = sin(a/2) 



(117) 



For a ^ 1, we therefore have ipc ~ ds/dt. 

Choosing w = t + r* as our time parameter, we can 
foliate the spacetime into three-dimensional slices. Using 
this slicing, we can derive the invariant distance per unit 
time using Equation (112): 



ds/dv — ^/Ji 



^v/A 



vr 



aV/A 



(118) 
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Likewise, solving (104) for 9{v), holding 9{v) = 0, we 
obtain a null generator k^ = (1, 0(/u), 9, 0) at the caustic, 
where 9 = y/]Je^^^/\/4:{l — e"/^). Using the metric at the 
horizon, the deviation between generators is then found 
to be: 



V'c 



aV/i 



/^ 



VT 



yu/4 



(119) 



Both of these properties, derived independently from 
each other, point to a cone angle of 



TT — a, a = 2y/]i 



Vl - e^/* 



(120) 



and are in agreement with each other. For small fj,, the 
caustic angle is very nearly equal to n for all reasonable 
times |z;| > /i far from the small hole's event horizon 
(where the perturbation theory is valid) . Only very near 
the merger - i.e. in the regime |f | < /i that the perturba- 
tion theory cannot probe - does the caustic angle deviate 
significantly from tt. 

Of course, these results should be taken with a grain 
of salt. While the caustic angle does not depend on the 
spatial coordinates we choose, it does depend on the time 
slicing. By a clever choice of coordinates, one may make 
a take any value one desires. However, relations (115) 



and (117) hold irrespective of the coordinate choice, al- 



lowing the caustic angle to be related to the invariant 
distance per unit time and the deviation between genera- 
tors in any time slicing. Moreover, the integrated caustic 
length Lc = 8y/JI is a slicing-invariant quantity as well. 

It is worth emphasizing that the results derived in this 
section are independent of the particulars of the plunge. 
We therefore say that the caustic structure is universal - 
it depends on the mass-energy of the infalling black hole, 
and not on anything else. 

Before we end this section, a few consistency checks 
are in order. First the perturbation theory breaks down 
when the field of the point mass is strong - i.e. near the 
small black hole's event horizon. Since the size of the 
event horizon is proportional to /i, our results are only 



valid for 9^ :^ 9eh 



fi. Note, however, that the caustic 
'^/^ so that in the extreme mass- 



horizon scales as 9c ^ fJ." 
ratio limit (/x ^ 1), 9c ^ 9 eh- Put in words, for a 
small infalling mass, the event horizon is tiny compared 
to the caustic horizon, and the field of the infalling mass 
is weak for rays of 9+ ~ 9c- Most caustic-forming rays, 
therefore, can be accurately described using the methods 
of this section. 

Second, in order for geodesic perturbation theory to 
be valid, we must require that (5r <C 1. At first glance, 
this appears to fail by virtue of the logarithmic term in 
(103 1. However, a closer inspection reveals that in the 



limit of small fj,, this term is of order unity only for rays 
of 9+ = 0(e"^/^) < 0{fi). As explained above, the 
gravitational field of m is strong for 9+ — 0{n), and 
the theory breaks down anyway. In the region where the 
theory does apply (0+ 3> /i), the radial perturbation is 
small everywhere. 
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Before Merger 
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FIG. 7: As black holes merge, the surface area of the large 
hole's horizon increases. This is due to rays entering the hori- 
zon through the caustic (AAc, shown in red), expnsion of rays 
near the caustic (Aj4„, blue), and expansion in the bulk of 
the horizon {AAb, gray). 



V. BLACK HOLE AREA INCREASE 

The surface area of a black hole of mass M is IGttM^. 
For an infinitesimal mass increase, the surface area in- 
creases by AA = 327rMAM. Thus, for the case in ques- 
tion - a point mass fj, <^ I falling into a black hole of 
mass M = 1, we expect the black hole's event horizon to 
grow by AA = 327r/j,. Naturally, we are inclined to ask: 
how much of this area increase is due to rays which enter 
the horizon through the caustic, and how much of it is 
due to the expansion of the horizon itself? 



A. General Principles 

Without loss of generality, assume the small black hole 
reaches the horizon at the +z axis - that is, with 9^0. 
We can then write the horizon area increase as a sum of 
three parts: A^l = AAc+AA^+AAb. Here AAc refers to 
the area increase due to rays entering through the caustic 
and AAn refers to the expansion of area elements in the 
"neighborhood" of the caustic - i.e. for ^ 1. The last 
term, AA^, refers to the increase due to the expansion of 
the bulk of the horizon. Figure [7] illustrates our point. 

A number of general principles allow one to solve the 
area increase problem with minimal effort. First, as men- 
tioned above, the total area increase must be 327r/i. Sec- 
ond, infinitesimal area elements can expand, but they 
cannot contract. 1 10] Third, area elements in the bulk of 
the horizon are minimally affected by the point mass's 
gravitational field, and therefore neither expand nor con- 
tract. 



B. Detailed Calculation 

The purpose of this subsection is to verify the points 
made above with explicit calculations of the three area 
contributions. 
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1. Area due to generators entering through the caustic 

Recalling that rays for which 0+ < Oc — Zy^ enter 
through the caustic, the area increase due to the caustic 
is simply equal to the area spanned by these generators 
on the final horizon, which is equal to 

AAc = 7r(26lc)2 ^ i67r^ (121) 

This is half of the total area increase. 

2. Area increase in caustic vicinity 

Consider, in the neighborhood of the caustic, an an- 
nulus 9 € [0+,0+ + d0+] at future infinity {v = +oo). 
Tracing these rays back to past infinity, we end up with 
a new annulus defined by: 



[9^,9^ +d0^] 
[9+-S9id+),t 



d04 



(122) 
-d0+)]{122,) 



where 69 {9+) is the 9 defiection of the ray due to the small 
hole's gravitational field. Assuming that the annulus is 
thin, i.e. d9^ <^ 9^, we find the new annulus has radius 
and width: 



d9. 



d9. 



859 
~d9 



(124) 
(125) 



At past infinity, the area of the annulus is A_ ~ 
8TTsm{9-)d9-. After the merger, the area grows to 
A_|_ = 2(2 + Sr)^TT sm{0 ^)d9-^-. The infinitesimal area 
increase is given by: 

6 A = A+-A^ 

= 8n [sm{9+)d9+ - sm{9^)d9^ + 6rsm{9+)d9+] 

(126) 

This is a general formula that applies equally to the caus- 
tic neighborhood and the bulk. 

Specializing to the caustic neighborhood, we set 9^ <C 
1. Equation (126) describes the area increase as due to 



two components - radial deformation, and angular de- 
formation. The radial deformation is trivial - at future 
infinity, the radius of the black hole increases by 6r — 2/i. 
In Section [TV] we calculated the angular deformation to 
be: 59[9^) = 4^/0+. Applying these substitutions, we 
arrive at the following area increase: 



5A = [l287r/^2gi 




(127) 



radial 



The angular part dies off quickly with increasing 9 and 
therefore only contributes in the neighborhood of the 
caustic. Ostensibly, it appears to be a second-order effect 



in to; however, when we integrate it over the neighbor- 
hood, we find an total area increase of: 

"'"'l287r^2^^167r/i (128) 

where 9c ^ ^fmax ^ 1 so that we only integrate over the 
neighborhood of the caustic. 

The neighborhood of the caus tic h as an area A^ ^ 
^max ^ 1- The radial term in (127) will increase the 
area by an amount 0{^An) <C 0{fi), which is much less 
than the area increase computed in (128). Therefore, 



to leading order in m, the area increase in the caustic 
neighborhood is: 

AAn = IG-TTfl (129) 

This contributes the other half to the total increase 327r/z. 

3. Area increase in the bulk 

Next, we turn to the bulk - that is, generators with 9 > 
^maxj with ^inax defined as above. Like near the caustic, 
an area element in the bulk will expand or contract due to 
two factors: angular deformation and radial deformation. 
The extent to which it does so is governed by (126). In 



the bulk, however, rays are weakly deflected, i.e. 69 <C 



f?+. This allows us to simplify ( 126 ) to the following form 



M = A 



cot( 



869 
~d9 



(130) 



(Here A is the size of an infinitesimal area element. It 
does not matter whether we use A^ or A_, since in the 
bulk of the horizon, A+ = A^ up to a relative error 

In the bulk, furthermore, the odd-parity perturbations 
cannot be neglected as they may in the caustic vicinity. 
Nevertheless, their effect on the horizon area is imma- 
terial. A quick glance at Equations ( 55][57 ) reveals that 
the forcing terms for the odd-parity perturbations are 
solenoidal in nature - they can shear area elements on the 
horizon, but they cannot expand them. The even-parity 
perturbations give rise to expansion, while the odd-parity 
terms give rise to shear. It is the expansion of area ele- 
ments that chiefly concerns us here. 

As above, we set 6r = 2/i; hence, 6r has no angular de- 
pendence. In Part | V A\ we asserted that the infinitesimal 
area elements in the bulk neither expand nor contract - 
i.e. that 6A = everywhere in the bulk. This gives a 
first-order ODE for 69{9). Solving subject to a continuity 
constraint at 9 — tt, we expect an angular deformation 
of the form: 

69 = 2^i^t_^ + (odd terms) (131) 

We can obtain precisely the same result by invoking the 
impulse approximation, applying the forcing terms in 
Equations ( [95|[96| to the bulk, and adding the odd-parity 
perturbations (which do not affect the area increase) . 
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VI. CAUSTIC FORMATION WITH STRINGS 

In the preceding sections, we have constrained our 
analysis to the infaU of point particles in the extreme 
mass-ratio case. For point particles, we found that the 
event horizon develops a caustic and that the caustic ac- 
counts for 50% of the black hole area increase. On the 
other hand, it should be eminently clear that extended 
objects larger than the caustic horizon will not form caus- 
tics when they plunge into the black hole. Naturally, one 
is inclined to ask: what about objects extended in one 
or two dimensions? In this section, we study the infall 
of such objects and conclude that, under certain circum- 
stances, they do form caustics, but that such caustics 
account for 100%, not 50%, of the black hole area in- 
crease. 




FIG. 8: A diagram of two "strings" falling into a black hole, 
one parallel to the horizon and the other orthogonal. In real- 
ity, most objects would fall in at oblique angles. 



A. Parallel Strings 

Recall our analysis for an infalling point particle in 
Section [TV] In the vicinity of the caustic horizon, 9^1 
and the event horizon's geometry is essentially Euclidean. 
We may therefore write 9 and (/) in terms of Euclidean 
coordinates x'' — {xi^x^) and recast Eqs. ( 97p9 l as: 



where 



d 

dv 



1 



6r ^ i^'■(x^w) + C(°)(^')-2(i2+i2) 
-2{x\ + xl) (132) 



d_ 
dv 



dv 



F''(x\v) = ~'^x'' 5{v) (133) 

x 



where x — \\x°'\\ and x is the normalized vector. 

Let us spread out the point mass into a string parallel 
to the event horizon. Since we are working to first order 
in the perturbation, our angular equations are linear, and 
thus obey the superposition principle. (The radial equa- 
tion is not linear, but once we have i", it is easily solved). 
We can therefore rewrite the forcing term in ( 133 ) as: 



F''{x\v) = 



q 



(134) 



where 5" = x"-—x°'' is a vector pointing from the mass ele- 
ment d/i to the observer at point 2;", and the vectors with 
hats {x°' , q"") represent normalized unit vectors. Solving 
the equations of motion prior to merger, we find that the 
displacement of any ray due to the string will be given 
by: 

5r = e''/4Ar + Cstr(w)- J(Aa;'^)2(e"/2-e"/4)(135) 



c^ + Ax''(e''/^-l) 



(136) 



Ar = - I A\og{q)d^i-AiJLC (137) 



Aa;° 



(138) 



and Cstt{v), defined in the same way as (1091, merely 



shifts the caustic in the r-direction and has no effect on 
its internal structure. 

Both Ar and Ax" obey the superposition principle, 
and can likewise be related to the structure of the string 
by analogues of Gauss's Law and Poisson's Equation: 



dad°-l^r = dal^x"- = 87rp 



(139) 



where p is the amount of mass which falls through the 
horizon per unit area in the (x^, x^) coordinates. The in- 
tegral J pdxidx2 would then correspond to the enclosed 
infalling mass. 

Now consider an arbitrary string falling in parallel to 
the horizon. Draw a Gaussian surface as per Figure [9] 
around any small line element of the string. A simple 
application of Gauss's Law allows one to relate the dis- 
placement vector Ax" on one side of the string to its 
value on the other side: 



Ax" 



SnAh" 



(140) 



The displacement quantities Ar and Ax" then take the 
following form near the string: 

Ar = A- Bax" - 47rA|x| (141) 

Ax" = S"+47rAsgn(x)n" (142) 

where A is the linear mass density of the string and A 
and B" are constants of integration. Note, critically, that 
unlike in the case of the point particle, the displacement 
field Ax" is roughly constant in the vicinity of the string; 
it does not depend on the distance between the string and 
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Ax! 




FIG. 9: A Gaussian surface chosen around a string falling into 
the event horizon. 



the observer. This displacement will cause the horizon 
generators near the string to converge towards it, forming 
a one-dimensional, linelike caustic. 

Consider again the Gaussian surface in Figure |9] The 
widths W-i- and W- on each side of the caustic are fixed 
so that Ar takes the same value at both edges of the 
Gaussian surface. Since they have the same Ar value, 
the rays at the edges of the surface will meet if and only 
if they cross in angular coordinates. This will happen if 
the total width w is less than SttA, the "jump" in Aa;° 
across the string. Therefore, the generators entering the 
horizon through the caustic form a small strip, containing 
the string, of width w = SttA. Integrating over the whole 
string, we find the total area increase due to the caustic: 



AAc = 4(87r^) = 327r/i 



(143) 



Thus, the linelike caustic of the string is responsible for 
100% of the increase in black hole area. This contrasts 
sharply with the 50% figure found for point particles. 
As we derived in Equation (142), the displacement field 
Aa;" is approximately a constant in the vicinity of the 
caustic, as opposed to the 1/6 dependence found for the 
point mass case. A constant displacement field makes it 
impossible for area elements near the caustic to expand 
(except at the edges of the string, which contribute neg- 
ligibly to the area). By our reasoning in Section^ none 
of the area increase comes from the bulk to first order in 
ji, so the entire 327r/i area increase must be due to rays 
entering the horizon through the caustic. 

Technically speaking, the line at which the null rays 
meet is not a caustic, but a crossover set — a set of points 
at which non-neighboring rays meet. The crossover set 
traverses a line inside the small strip comprising the caus- 
tic horizon, and is terminated on both ends by caustic 
points — points at which neighboring rays meet. This is 
to the description given by Husa and Winicour [TT] for 
asymmetric mergers. By extending the compact mass 
onto a stringlike object, we break the rotational symme- 
try of our problem, which in turn causes the horizon to 
be perturbed asymmetrically, forming caustics as well as 
a crossover set. 



B. Orthogonal Strings 

Consider instead a string which falls into the black hole 
orthogonal to the event horizon. In this case, returning 
to polar coordinates, we simply spread out the delta- 



function in (98), obtaining: 



d 

dv 



dv 



e 



fiv) 



(144) 



where f{v) integrates to unity. Now suppose that the 
string is very thin, so that the forcing term f{v) is very 
small. In this case, we may neglect 9{v) as small com- 
pared to 0(v), yielding the following equation of motion: 



difl 
dv 



Sfifiv) 



(145) 



Rays will cross 9 — and form caustics if and only if: 



9+< ^ 



(146) 



This is slightly larger than the caustic horizon found for 
a point particle. It corresponds to a black hole area in- 
crease of: 



AAc = 4(87r^) = 327r^ 



(147) 



Again, the caustic is responsible for 100% of the black 
hole area increase, just as in the case of the parallel string. 
Other objects, including oblique strings and two- 
dimensional sheets, may be studied using the methods 
of this section. Like parallel and orthogonal strings, such 
objects are expected to produce caustics which account 
for 100% of the black hole area increase. The analysis, 
however, is less enlightening and is not presented here. 



VII. CONCLUSIONS 

In this paper, we used black-hole perturbation the- 
ory to model the event horizon of an extreme-mass-ratio 
merger. While it is straightforward to obtain the event- 
horizon deformation by numerically integrating the Zer- 
illi and Regge- Wheeler equations — and then geodesic 
equations, we have taken advantage of the extreme mass 
ratio and used an impulse approximation, which allowed 
us to find a universal geometry of the event horizon for 
all such mergers. 

While most of the large hole's horizon generators orig- 
inate on the future horizon and remain on it as the black 
holes marge, a small subset originate on the past horizon 
of the large black hole, wander away, and are gravitation- 
ally lensed by the small black hole's field. They subse- 
quently enter the horizon at the caustic and remain on 
the horizon after the merger. The bundle of rays entering 
the horizon through the caustic is given by: 



< v/4m 



(148) 
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where, without loss of generality, we assume the small 
black hole reaches the horizon at the z-axis. 

We find a caustic which is qualitatively similar to that 
seen in previous studies of head-on collisions. Locally, 
the horizon is cone-shaped at the caustic, and the cone 
becomes increasingly pointed as the small black hole ap- 
proaches the horizon. We find that the cone angle can be 
related to the invariant length per unit time in a given 
time slicing, and the total invariant length is calculated to 
be i^/Jl■ We find that the structure of the caustic is uni- 
versal, depending only on the masses of the two merging 
black holes. While the trajectory of the small black hole 
may affect the deformations in the bulk of the horizon, it 
does not materially affect the caustic. On the other hand, 
cone-shaped caustics are not the only caustics that may 
be formed; if it were possible for a string-shaped object 
to fall into the black hole, it would produce two caustic 
points at the ends of the string, at which neighboring 
rays meet, and a crossover line at which non-neighboring 
rays meet. Extending the point mass into a string breaks 
the isotropy of the problem and necessarily changes the 
shape of the caustic. 

Lastly, we found that the black hole area increase is 
half due to rays entering the horizon through the caustic, 
and half due to the expansion of rays in the vicinity of the 
caustic. This property is also universal - while different 
infall trajectories may result in differing degrees of shear 
in the bulk, they all produce identical results at and near 
the caustic, where all of the area increase happens. 

That said, there are a number of limitations to this 
study. First, by restricting ourselves to first-order per- 
turbations, our results are only valid when the gravi- 
tational field of the small hole is weak. Thus, while 
weakly lensed generators are correctly described using 
this model, stronly-lensed generators - those with 6+ ^ fi 
- are not. Our perturbation theory cannot resolve the 
caustic at lengthscales comparable to the small hole's 
Schwarzschild radius, and it is plausible that at these 
lengthscales the conical structure of the caustic breaks 
down. In order to go beyond the results of this paper, we 
would need to either consider higher orders in the per- 
turbation expansion, or develop a qualitatively different 
technique for modeling the merger. 

Finally, while this paper is restricted to Schwarzschild 
black holes, most black holes in the real universe have 
spin. Moreover, there does not exist a straightforward 



formalism for metric perturbations in the Kerr metric, 
so the approach taken by this paper is not easily trans- 
ferrable to Kerr. Nevertheless, the dramatic simplifi- 
cations and the universal geometrical features shown in 
this paper indicate that, for extreme mass ratios and for 
near-horizon geometry, it may not be necessary to start 
with the full metric-perturbation equations. It is conceiv- 
able that a Rindler Approximation near the Kerr horizon 
would allow us to demonstrate the same caustic geome- 
try (in a coordinate system co-rotating with the horizon) 
as in a plunge into Schwarzschild black holes. 
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Appendix A: Even-Parity Perturbations 

In the following two Appendices, we give the perturbed 
Einstein Equations and discuss how the metric perturba- 
tion is related to the Zerilli and Rcggc- Wheeler functions. 
To do so, we use the formalism and results of Zerilli's pa- 
per [31] • 



1. Einstein Equations 

The perturbed Einstein Equations take the form: 

SG^, = SttT^, (A1) 



We express the energy-momentum tensor in terms of 

1(0) 4(1 



seven spherical harmonic components: Aim, ^i^, A, 



Bim, -B;„, Fim, and Gim- 

The tensor, expressed in matrix form, looks like the 
following: 



1 

V2 



v24:]y 






■i 4(-'-) v'™ 
1 -» 



V2A 



Im 



rCiBiraY'^ 

r^{Gi„,Y'"' + C2Fim 



W 



lni\ 



trG.Biy^^ 
rGiBimY'^ 

r^ Sin^ 0{GlmY''^ + C2FirnW'" 



(A2) 



Here we have defined 



Ci 



V^i\ + iy 



Co 



2VA(A + 1) 



(A3) 
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and 





' d ^ 
,90 y 


yrlm 


^Im ^ 


902 




1 9^1 
sin2 902_ 



F 



/m. 



(A4) 

(A5) 



The symbol * stands for terms obtainable through symmetry of T^i,. Seven distinct Einstein Equations are found in 
all: 



-8nA 



(0) 
Im 



{r-2yd^K (r-2)(3r-5)9X (r - 2)^ 9^2 r-2 



dr 



dr 



{H2 - K) 



(A + l)(r-2) 



Stti (1) ^ d_ 
72 '" dt 



dK K - H2 
dr 



r r(r — 2) 



K 



A + 1 



Hi 



{H2 + K) 
(A6) 

(A7) 



92 if r-l dK 2 (97Ji 1 9iJo 



8TTtV2rCiB" = 



(r-2)2 at2 r(r-2) ar r - 2 dt r dr r{r-2) 



(0) r-2 dHi , 2 ^^ a(i/2 + K) 



{H2 - K) 



A + 1 

r(r-2) 



(i^ - i^o) 



dr 



-^Hi — 



dt 



8.V-2Ciir-2)B,„. = -^ + — ^^^^ + 4^o + ^(i/. - i^o) 



8TrV2C2r^Fi 



dt r 

Ho — H2 



im 



/- _ r d^K r-2d^K 2{r - I) dK r d'^H2 d^Hi r-2d'^Ho 

™ r — 2 9f2 ^ 5^2 J.2 Q^f, r — 2 dt'^ drdt r dr"^ 

2{r-l)dHi r^l8H2 r + ldH„ A + 1,^ ^^ 
4 —[H2 - Ho) 



r{r - 2) dt 



dr 



J.2 Qj. 



(A8) 
(A9) 

(AlO) 

(All) 

(A12) 



2. The Zerilli Function 



Instead of demonstrating how the Zerilli function is motivated, we simply write down its expression, 

.2 ^_2 



Z 



-K 



Xr + 3 Ar + 3 



H,dt 



(A13) 



and outline how to show that it satisfies a wave equation, and how to obtain all perturbation fields from it. Through- 
out this section, we assume integrals in time range from —00 to t. We also assume all metric perturbation fields 
{Ho,Hi,H2,K) to vanish at —00. 

In order to arrive at the wave equation, we will have to take r derivatives and t derivatives. As it turns out, r 
derivatives can be simplified when we repeatedly apply 



dK _ r-3 

dr r{r — 2) 



^+1^2 + ^-^ I ^1^^ - 4^*^/2 / A\^ldt , 



which can be obtained by integrating (AT) in time from —00 to t, and 
d 



„ , Hidt = / -^dt = — ^ — - / Hidt —(H2 + K) ^, , 

dr J J dr r{r ~ 2) J r - 2^ ' r-2 '" ' 



8TTiV2r^Ci^^o) 



(A14) 



(A15) 



which can be obtained by integrating and simplifying (A9). Incidentally, after substituting (A14) and (A15), dZ/dr 
may also be expressed in terms of K and / Hidt: 



K + - ' v: '::: ■ - / H,dt ^f^^ / A\'jdt - --- -^ / B^yt . (ai6) 



dZ _ r(Ar2 - 3Ar - 3) A(A + l)r2 + 3Ar + 6 

Ih^ ^^ (r-2)(Ar + 3)2 ^ r(Ar + 3)2 

Further calculation shows that d^Z/dr^ is a combination of K, J Hidt, H2 and source terms 



Xr + 3 J '™"^^ " Ar + 3 
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On the other hand, d^Z/dt^ can be expressed also in terms of K , J Hidt and H2, with the help of Eq. (A8I, plus 



dH2 
dr 



r-3 



r - 4 



r dHi A + 1 



r(r-2) r(r-2) 



2 dt 



Hidt 



dir^Firn) r{r - 3) 



dr 



r-2 



Fi, 



8nV2CirBi, 



(A17) 



which can be obtained by substituting (All I and (A14) into (AlO) 
Putting together d^Z/drl and d^Z/dt"^, we obtain 



92 d^ z _ z 



Im 



(AI8) 



where 



Vi^ = (1-2/r 



2Y{\ + l)r^ + eAV^ + 18Ar + 18 
r3(Ar + 3)2 



2 _ 87r(r-2)2 87r(r - 2)^ 8^7r(r - 2) 8N/2^i(r - 2)(r(A + 3) - 3) 

"" " Ar + 3 ^'™+yA+T(Ar + 3) '" v/A(ATI)"" r(Ar + 3)2 

87r»(r-2)(r2A2 + 3r(A-2) + 12) f („) 4727rz(r - 2f f (1) 

+ VA + lr(Ar + 3)2 i "^'"^ '^ " A. + 3 7 ^"" -^' 



87ri(r-2)2 
\/ATT(Ar + 3y 



5i:iVc;^ 



(A19) 



4'2dt 



(A20) 



Once the Zerilli function is known, the metric perturbations may be calculated accordingly. Inverting Equations ( A13) 
and (AI6I, we find: 



K 



H, = 



r2A(A + l)+3rA + 6 r-2dZ 4:V2nir{r - 2) /" .(1) , 8%/27rir(r - 2)C: 



r2(Ar + 3) 



-Z 



r dr 



Ar + 3 



A 'dt + 



3 + rA 



'- / Biyt (A21) 



^A-3rA-3aZ , d'^Z , 4V2TTir\(^i) 8V2TriCir\(^Q) 



(r - 2){\r + 3) dt drdt Ar + 3 



-A 



Irn 



rA + 3 '" 



(A22) 



Substituting these formulae for Hi and K into ( A14), we find an expression for H2 in terms of the Zerilli function: 



Ho = - 



^A2(A + l) + 3r2A2 + 9rA + 9^ r^A-rA + S^Z 



d^z 



r2(Ar + 3)2 

8V2mCir{r^X{X - 1) + 6r(A - 1) + 15) 
(Ar + 3)2 



r(Ar + 3) dr dr^ 



B\yt 



Ar + 3 / '™ ^'■ 



4V27^^r(r2A + 6r - 6) /" .(1) ,^ , 4y27rzr2(r - 2) f (1) 

(Ar + 3)2 J ^"»^^+ aTTs y ^'"^-''^ 



(A23) 



3. Source Terms 



which determines the geodesic motion via 



The point mass of the small black hole traces out a 
trajectory (T(r), i?(r), 6(r), $(t)) in the Schwarzschild 
spacetime. Because of the spherical symmetry of the 
problem, we may set 6(r) = 7r/2 by an appopriate ro- 
tation. The trajectory of the particle is characterized by 
two invariants: the specific energy E = (1 — 2/r)M*, and 
the specific angular momentum L = r'^u'^ . The effective 
potential of the black hole is 



U{r)^{l-2/r){l + L'/r') 



(A24) 



{dr/drf + U{r) = E'^ 



(A25) 



For a rest mass of toq, the energy-momentum tensor 
of the particle is. 



T,. = ^f ^^'5(r - i?(t))j(2)(^ _ n(t)) (A26) 
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The even-parity content of this tensor corresponds to: 

Aim 



j^^^^E^ U[r]5{t ~ T{r))Y:^{t) (A27) 



^ - ^^^Sit-Tir))Y,Ut) 



Im 



ir 



(r-2) 



iV2CimmoEL 
Bim = -, ;;^ 6{t~T{r))Yi„^{t) 



r- 



•ir-2) 



b\'2 = 



V2CimmoEL{r - 2) 
r4^S2 -U[r] 

■moC2L^ 



Im 



V2r^y/E^ - U[r] 



(A28) 

(A29) 

S{t-T{r))Y;^{t) (A30) 

6it-T{r))WiUt) (A31) 



The spherical harmonic Yi*^{t) refers to the value of Yj^ 
at (0(i), <I>(i)), and is therefore a time-dependent quan- 
tity. The same holds for Wi*^{t). 



Appendix B: Odd-Parity Perturbations 

The odd-parity perturbations will give rise to three dis- 
tinct Einstein Equations, one of which may be discarded 
as redundant after the Regge- Wheeler gauge is imposed. 
There are two metric perturbations, ho and hi, which 
may be written in terms of a Zerilli-type function (al- 
though its existence predates Zerilli's work|29j) Q, which 
in turn satisfies a sourced wave equation. 



1. Einstein Equations 



As in the even-parity case, the perturbed Einstein 
Equations take the form SGf^^ — SttT^^. This time, how- 
ever, the energy-momentum tensor takes the following 
form: 



T, 



fiiy 



1 

71 



liO) 



Circsc(6')F™g 

* iCircsc{0)Y'^Qi, 

* * iC2r^cac{9)X^'^D, 



Im 



* * 



* 



-Cirsm(e)Y'j^Q\°2 
-iCirsmi9)Y'g"'Qim 

-iCW sin(0)X'" Am 



(Bl) 



As in the even-parity case, we have replaced some of the quantities with asterisks to simplify the expression. The two 
equations we choose to solve are the SGrtp and SGetp equations, which take the following form: 



OTT i^l^lm = - 



Stt- 



V2 
V2 '" " 2 



2(r-2) 



d'^hi d^ho ^2dho ^ 2A(r - 2) 
"9^2"" '[fi^t^ ^ 



r-2dhi 
r dr 



M- 



r dt 
r dho 



hi 



-2 dt 



(B2) 
(B3) 



2. The Regge- Wheeler Function 

The Regge- Wheeler function can be expressed in terms 
of hi as: 



r-2 
Q=^^hi 



(B4) 



Like the Zerilli function, Q obeys a simple wave equation 
with an analytic potential and source term. The equation 
is given by: 



UQ-V^^Q^S^^^ 



where 



u'- 



1- 



2(A + 1) 6 



(B5) 
(B6) 



S?ra = -8V27rz(r - 2)C2 



d fr- 
dr \ r 


^D, 


V2A(r 


-2) 



Q 



Im 



(B7) 



By solving this wave equation, we find the Zerilli func- 
tion, but what we are really after is the metric perturba- 
tion. Substituting (B4| into ( B2||B3 ), we may relate hi 
and ho to Q: 



ho = "^ fQdt+{r-2) f^dt 
+ 8iV2nr{r - 2) / Dimdt 



r-2 



Q 



Source Terms 



(B8) 
(B9) 



Following the methods of Section |A3| we can obtain 
the source terms for the odd perturbations. It is only 
necessary to obtain Dim and Qim, which are provided 
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FIG. 10: Plot of \Z{(jj)\ versus cu. Larger values of / are de- 
noted by darker lines. 



FIG. 11: Plot of \Ke(Z (ijj))\ versus oj. Larger values of I are 
denoted by darker lines. 



below: 



to the Zerilli-Moncrief function: 



D,,. 



^Im — 



\/2C2mmoi^ 



r^y/E^ - U{r) 
r2(r-2) 



(5(i-T(r))Y:^''' (BIO) 
5(i-T(r))y^"* (Bll) 



Appendix C: Behavior of Zerilli and Regge- Wheeler 
Functions 



The ZeriUi and Regge- Wheeler functions are com- 
puted by Fourier transforming Equations ( A13) and (B4) 
to the frequency domain and computing Z(r^,uj) and 
Q{r^, , ut) using a standard ODE solver like Mathematica's 
NDSolve [] . As per the standard practice, we impose an 
outgoing boundary condition at r* = oo and an ingoing 
boundary condition at r* = — oo. Since we are interested 
in the values of Z and Q at the horizon, we take Z^r^,) 
(or Q{r^,)) with r* sufficiently negative (i.e. r* — —20) 
that, for all practical purposes, it is on the horizon. Once 
we have the horizon Z{uj) (or Q{uj)), we Fourier trans- 
form it into the time domain, obtaining Z(v) (or Q{v)). 
Our Fourier transform algorithm takes into account the 
asymptotic behavior of Z and Q as w — )■ and oo. 



1. Zerilli EYinction 

In this section, we consider only the case of a radial 
infall. The Zerilli function for a general infall will be 
qualitatively similar. 

Before calculating Z{oj), let us try to determine its 
general structure using intuition and some asymptotics. 
Consider first, the case to ^ 0. This corresponds to the 
behavior of Z on long timescales. We know that for v ^■ 
— oo, Z ^ 0. To get the limit as w — >■ -t-oo, we relate Z 



Zm — 



A + 1 



■^ r-2 / ^ dK 

Ar + 3 V or 



(CI) 



The Moncrief function is a linear combination of the 
metric perturbations and, by the No-Hair Theorem, must 
vanish as w — > +oo. It is related to Z by: 



ij\/ 



4\/27r?>2(r-2) 
(A + l)(Ar + 3) 

(A-f l)(Ar + 3) 



Al'^dt 



(C2) 



Near the horizon, therefore, the Zerilli function takes the 
following V — > +0O limit: 



Z 



SnfirY* 



(A + l)(2A + 3) 



(C3) 



Thus, on long timescales, Z{v) is approximated by the 
Heaviside function Z{v) « Zrx,H{v). Taking the Fourier 
transform, we find the lower asymptotic limit: 






(C4) 



Now consider the asymptotic behavior as w — >■ oo. This 
encodes the short-lengthscale behavior of the Zerilli func- 
tion near the path of the particle. On these lengthscales, 
we can assume that Z is sourced solely by a delta-function 
component along the particle's trajectory; this gives rise 
to a discontinuity in the Zerilli Function's derivative. 
Near the horizon, the source term for Z may be writ- 
ten in the following form: 



2A + 3 



6^^^[x'^-x'^{r)]dTY,l 



(C5) 



where x°'{t) refers to the path of the particle. Since 
the potential V{r) vanishes near the horizon, the Zerilli 
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Equation reduces to a free, sourced one-dimensional wave 
equation. Solving this by the method of characteristics, 
we find the following jump in the derivative dZ/dv: 



AZ' 



dZ, 



Im 



dv 



dZ, 



Ira 



dv 



"2A + 3 '" 



(C6) 



This gives rise to the following asymptotic behavior: 



1 ^TimE^ 
/2^a;2 2A + 3 ^ 



:y: 



m 



as u) 



(C7) 



The frequency-domain Zerilli function may be computed 
numerically by solving the equation: 



drl 



Z + uj^Z -V^Z = S^{uj) 



(C8) 



We plot this function in Figures [TO] and [TT] and compare 
it to the asymptotic limits derived above. 

We wish to derive a scaling relation for the Zerilli func- 
tion for large I. The first step in doing so is to note that 
the jump in Z'{v) scales as l^'^Yi*^ , whereas the long- 

timescale jump Z^o scales as ?~^1^;„ • For large I values, 
it is the contribution of the former that dominates. This 
contribution may be extracted from Z(uj) by taking the 
real part. This is plotted in Figure [TT] 

The limiting value Zq = Re(Z(w — 0)) may be ob- 
tained by looking at the asymptotic behavior of (C8| 
Since F - P ^^^ j^gj-^) 

large I, Zq ~ 1-%T 
lines in Figure [TT] 

The cutoff frequency for |Re(Z(a;))| may be obtained 
by equating the asymptotics, which scale as /"''i^;^ and 
^-2^-2yjn respectively. We find ujc ^ I, or correspond- 
ingly in the time domain, tc ^ l~^. Taking the Fourier 
transform of Z{uj), the value of Z{v — 0) is proportional 
to the area under the curve of Z{uj). This area, in turn. 



l-^Y, 



.(0) 



Im 



it follows that for 



This behavior is plotted in dashed 



scales as Zqujc 



1''y;^'\ so Z{v = 0) 



l-^Y, 



.(0) 



We 



thus obtain a scaling relation for Z(v) for large I: 

Zlm 



,_3^*(0) 



fiv/l-') 



(C9) 



where /(^) is an undefined function which is independent 
of V. Plotting Zim/l~^Y,*J; against v/l~^, we obtain 



the graph in Figure 12 Empirically, Z{v) approaches an 
exponential: 



Zu 



Sirr 



iV*W„~l\v\/2 



(CIO) 



What effect does this Z{v) have on the horizon? Recall 
that the metric perturbations act on the horizon through 
the forcing terms ( 55][57 ). For even perturbations, all 
forcing terms are proportional to 
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(e) _ 



32TnnoEY 
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Im 
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FIG. 12: Plot of Z{v)/l-%*J^°'' versus v/T^. Larger values 
of I are denoted by darker lines. The limit curve is shown in 
red. 



Note that /;^ contains two delta-functions - one explic- 
i tly, a nd one due to d^Z/dv^. Noting the asymptotic form 
( CIO I for Z, we can see that both of these delta-functions 



precisely cancel out! The effect of the Zerilli function is 
to smooth out the delta-function in /^ ^ , resulting in a 
forcing term which is continuous and has a characteristic 
timescale tc ^ l~^. If we look on timescales much longer 
than l~^, the forcing term may be approximated by a 
delta- function, as we showed in Section IIIIl 



2. Regge- Wheeler Function 

Like the Zerilli function, the Regge- Wheeler function 
Q satisfies a wave equation with a very similar poten- 
tial. The source term (B7), however, is slightly different 



- unlike the Zerilli function source, the Q source does 
not contain a Heaviside-function component, and it does 
contain a 5'{t — T{r)) component. These differences will 
make the asymptotic behavior of Q{uj) slightly different 
from that of Z. 

Like Z{uj), Q is computed using a Mathematica ODE 
solver script. We plot Q, rms averaged over m values, in 
Figure [13] 

We proceed to determine the general asymptotic be- 
havior of Q. It will turn out that Q's contribution to 
the forcing terms is negligible near the caustic, so we do 
not take the time to derive exact formulas. For large 
I, the general behavior of Q may be determined by an 
asymptotic analysis of the Zerilli Equation ( B5 ) : 



dr'i 



Q + {lj^-V^)Q = SQ 



(C12) 



In the limit oj — > 0, the dominant term in S^ scales 



I '^Y,l„ a, while V'^ ^ P. Since these terms 



as S'^ 

dominate in the limit of large I, we find 



Q' 
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r'^Y, 



Imfi 



as uj 



(C13) 
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FIG. 13: Plot ol {Q{oj)) versus u. Larger values ol I are 
denoted by darker lines. 



Likewise, for large oj, we must balance uj'^Q on the left- 
hand side against S*^ on the right-hand side. The dom- 
inant part of S'^ in this case scales as l~'^mujYi*^ g ^ 
l^'^ujYi*^^ „. Thus, the appropriate scaling for Q is: 



hyifi 



-Y, 



Imfi 



as w ■ 



(C14) 



Like in the even case, the cutoff frequency for Zodd{^) 
scales as Wc ~ ^ - or conversely in the time domain, the 
odd-parity Zerilli function has a characteristic timescale 

The odd-parity forcing terms from ( 55]|57 ) are all pro- 
portional to: 



Ao) ^ 

■> Im 



dv 



1 



(o) 



Q 



(C15) 



Like its even-parity equivalent, f^^ has a characteristic 
timescale of ic ~ l~^- If we look on timescales much 
longer than l~^, the forcing term may be approximated 
by the delta-function derived in Section |III[ The odd- 
parity term turns out to be sub-dominant, and can ac- 
cordingly be neglected. 



Appendix D: Validity of Approximations 

1. Delta- Function Approximation 

Recall from Appendix |C 1| that the Zerilli function 
scales as a function of v/l^'-, and therefore for large I, 
the forcing terms increasingly resemble delta-functions 
in time. The impulse approximation, which treats these 
terms as delta-functions, results in a deflection of order 
A9 ~ ^^/^ for rays near the caustic horizon; as we show 
in this section, relative error in 9{v) incurred by imposing 
the impulse approximation is, at largest, of order ^^''^. 
This is only an upper bound, however, and empirically 



the relative error in 9{v) appears to be of order /x^/^. 
Therefore, in the extreme mass-ratio case /x <C 1, devia- 
tions from the approximation may be safely ignored. 



a. Proof by Asymptotics 



cie) 



First, we neglect the forcing term components /j^ with 
I < /i^/^ and show that this introduces an error of order 
/i. Next, considering only components with I > fJ.^''^, we 
show that the forcing terms may be written in terms of 
a smoothed-out delta function f{v) with timescale r < 
/i^''*. Deviations from the delta- function approximation 
give rise to an error at fx^^^r < ii^^^. 

Consider a particle falling in at the z-axis. Recall from 
Section|TTl]that for a given /, the even-parity forcing term 
component /,^ takes the form: 



— &'>Y 

TTmoEY*„^ 
A + 1 



(0) 



-Yim.eSi{v) 



(Dl) 



where 6i is a delta-function spread out over a time in- 
terval of order t; ^ l^^ which integrates to unity. First, 
we show that for generators in the vicinity of the caustic 
horizon 6 ~ H^^'^, it is possible to ignore all terms with 
I < /i-i/4 Re-expressing Yio{9,(j)) - ^^{21 + l)/4T:Joil9) 
for 6* ^ 1 (up to a relative error 0{9) ~ 0{fj}^-^)), we can 
obtain the following expression for the forcing term: 



Fff 



E 
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^{21 
21(1- 

H{21 
2(/ + l) 



i^Mium 
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5i{v)J'^{l9) 
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Now for / < 1/6* - ^-1/2^ J'^{19) - 19/2 - m}^' I, giving 

3/2 

a contribution to Fg of mg I. Summing over all spherical 

harmonic terms with / < m^ , we find a total contri- 
bution to Ffl of: 



F, 



i<j. 



M 



3/2 



^^ 



-1/4 
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(D3) 



The total forcing term Fq as derived in (92 1, on the other 



han d, is of order ^/9 • 



D3 



{I < M"'/': 



Therefore, the contribution 
terms is negligibly small. We 



of "small-r' 

may therefore assume, to leading order in /i, that only 
terms with I > ^~^/^ contribute to the structure of the 
caustic. 

The forcing term, as a function of f , is thus concen- 
trated within a small time interval of order r < ^^'^, 
and may be written to leading order in m as Fg — 
— {^/9)f{v), where f{v) integrates to unity. Write f{v) 
as a delta function plus a residual term: 6{v) + ^f(v), 
where e = 1 is a placeholder. Then 9{v) may be ex- 
panded in terms of e: 9 = 9^ + e9i + . . .. Geodesic equa- 
tion (Esl, written as 9" -9' /A = -{^/9){6{v) + efiv)) for 
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the current forcing term, then breaks down into a series 
of equations: 



%--0', = -^S{v) 



(D4) 



^i-7^i = ^O,6{v) + F{0o)f{v) (D5) 

4 «o 

We are not interested in the exact forms of 6*0 and 9i. 
Rather, we wish to determine how they vary to order of 
magnitude. In the time interval r < ij}' over which the 
ray is lensed, we find that: 






fJ- 



1/2 



r^^i/V 
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Likewise, 9i behaves as: 




FIG. 14: Plot of the horizon generators as lensed by Eq. |D9| 
for three masses: /^ = 0.1 (light gray), 0.01 (dark gray), and 
0.001 (black), and the delta-function (red, dashed). 
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Thus, 9'i ^ fi^'^T. This is a factor r < fi^'^ smaller than 
the deflection obtained using the delta-function approx- 
imation. It follows that all deviations from the impulse 
appoximation are smaller than the delta-function result 
by a factor of r < /^^''*, and therefore the impulse ap- 
proximation is valid for small /i. 

The same does not hold for the radial forcing term. 
As we show in equation [95) the radial term F^ scales as 
O(TOlogm), while the low-l contributions to this term 
scale as 0(rn). However, the \aw-l contributions are 
constant on the 0{rn}^'^) lengthscale associated with the 
caustic; this is best illustrated by noting that the spatial 
derivatives of F"^ are proportional to F^ and F"^, which, 
as we shoed in this section, are well appoximated as im- 
pulses. Thus, the impulse approximation for F^ is valid 
up to a constant term C^^' [t) , which depends on time 
but not position, near the caustic. This constant term is 
a function of the low-l perturbations and is 0{m). 



b. Empirical Verification 



We can verify this empirically using the shape of the 



impulse, as derived in Eq. (84) of Section IIIC 



F'{e,t) 



'4(6l2+t2/4)3/2 
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To assess the validity of the impulse approximation, we 
propagate generators using two different forcing terms: 
the delta- function in (87) and the spread-out impulse in 
(D9). This is illustrated in Figure 14 where we plot the 
paths of the generators for three masses: ^ = 0.1, 0.01, 



and 0.001. As the particle's mass becomes smaller, the 
delta-function approximation becomes more accurate. 

For the two impulse shapes, we compute the displace- 
ment 59 = 9\_^ and compare the two displacements. 
Defining a relative error 



\59i 



we find empirically that this error scales as: 



(DIO) 



(Dll) 



Therefore, for rays near the caustic horizon, the relative 
error accrued by imposing the impulse approximation is 
0(^^/2)^ -^^ell below the upper limit of fi^^^ derived in the 
previous section. 



2. Neglecting Odd Perturbations 

In general, the horizon generators are deflected by two 
types of forcing terms - odd-parity and even-parity. For 
rays near the caustic horizon, we show in this section that 
the ratio of odd-parity forcing to even-parity forcing is at 
most 0(/^^/^). Thus, for large mass ratios, the odd-parity 
terms may be neglected. 

The forcing terms may be written in terms of /^^ and 
Ao). 
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Accordingly, for a given Z, the relative importance of odd 
and even perturbations is given by: 



Qdd;^ 
Even;™ 
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Jim 
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(D14) 
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Therefore, /;^ < //^ for all I. Consider a null generator 
in the vicinity of the caustic horizon - 6* ~ /i'^/^. As 
we showed in the previous subsection, the even-parity 
solution is dominated by terms with I > fi~^'*. The 
even-parity contribution for I < /i^^'^ is of order /i^" 
smaller than the total contribution, and therefore can be 
neglected. Since /;^ < //^ , it follows that the odd-parity 
contribution of terms with I < ii~^'^ can be likewise 
neglected. Odd terms for which I > /i^^/* can likewise 
be neglected to leading order in fi because fi^/fl^ '^ 



I 



"'<m'/^- 



It follows that the odd perturbations as a whole may 
be neglected. The total odd-parity contribution to the 
forcing term is insignificant compared to the even-parity 
contribution: 



Odd 

Even 



< m'/' « 1 



(D15) 
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